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SUMMARY 

The coupling between vibrational nonequilibrium and dissociation has been shown 
previously to have a large effect on flow-field properties when dissociation is significant. 
However, most of the flow-field analyses presented to date used measured dissociation- 
rate constants in which the vibrational energy was assumed to be in equilibrium where the 
measurements were performed. Treanor has shown that this is not the case and that 
appreciable errors may result from the assumption of vibrational equilibrium in the 
measurements of dissociation- rate constants. 

To illustrate the effects of vibrational relaxation on the measurements of 
dissociation- rate constants, this report presents computed flow-field properties for the 
nonequilibrium flow of oxygen over a wedge using both the coupled vibration-dissociation 
and the coupled vibration-dissociation-vibration models for two "limiting" conditions. 

One condition or interpretation assumes vibrational equilibrium in the rate constant 
measurements, whereas the other assumes that the rate constant measurements already 
include the vibrational nonequilibrium effects. The results showed significant differences 
between the flow-field properties for these two limiting conditions at a Mach number of 
27.17 (where dissociation is appreciable) but small differences at a Mach number of 9.45 
(where dissociation is small). 


INTRODUCTION 

The current emphasis on space travel has created an ever increasing demand for 
more accurate computations of the flow field surrounding reentry vehicles traveling at 
altitudes and velocities where real gas effects are important. In the most general sense 
"real-gas effects" means the consideration of finite rate gas chemistry and the effect of 
coupling between the various energy modes on these rates. Even at moderate hypersonic 
speeds, where dissociation is significant but ionization negligible, it has been shown that 



the coupling of vibrational nonequilibrium and dissociation is important (refs. 1 and 2) 
even though the coupling mechanism itself is not very precisely understood. 

A number of methods have been suggested for calculating the effect on reaction 
rates of interactions between the vibrational energy states and dissociation. Hammerling, 
Teare, and Kivel (ref. 1) developed a coupled vibration-dissociation (CVD) model which 
accounts for the vibrational nonequilibrium effects on the rate of dissociation. This CVD 
model was modified by Treanor and Marrone (ref. 2) to account for the dissociation non- 
equilibrium effects on the vibrational energy as well as the effect of vibrational nonequi- 
librium on dissociation (CVDV model). Reference 3 presents computed results for the 
nonequilibrium flow of nitrogen over a wedge, where the CVD and CVDV models are com- 
pared with the uncoupled model. Additional refinements to the CVDV model have been 
suggested by references 4, 5, and 6; however, their application to flows other than one- 
dimensional problems is quite tedious. 

In most of the applications to date of the CVD and CVDV gas models, the computa- 
tions were performed by assuming that the dissociation- rate constants were experimen- 
tally measured in a region where the vibrational energy was in equilibrium with the 
translational mode. Even the experimental investigators of references 7, 8, and 9 
assumed that the vibrational energy was in equilibrium at the position aft of the shock 
wave where dissociation- rate constants were measured. However, Treanor, in a more 
recent report (ref. 10), indicates that the vibrational energy was not in equilibrium where 
the measurements were taken and that appreciable errors may result from the assumption 
of vibrational equilibrium. Also, Wray (ref. 11) found experimental evidence of the effect 
of vibrational relaxation on dissociation rates. 

This report presents the inviscid flow-field properties for the nonequilibrium flow 
of oxygen over a wedge using both the CVD and CVDV coupling models. In order to 
assess the importance of the vibrational relaxation effects on measurements of 
dissociation-rate constants, flow-field properties are computed for two different interpre- 
tations of the rate constant measurements, referred to as "limiting" conditions herein. 

One condition considers the rate constant measured with vibrational equilibrium (as 
assumed by other investigators), whereas the other condition assumes that the measured 
rate constant already includes the nonequilibrium vibrational effects. For any given 
coupling model, calculations using rates determined in any other fashion should lie 
between these two limiting conditions. It should be noted, however, that the two limiting 
conditions used here may not be proper bounds for all types of flow fields. 
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SYMBOLS 


a m>b m ,d m 

E(T,T V ) 

E v 


vectors used in equation (18) that correspond to terms in equations (1) to (5) 
term defined by equation (10) 
vibrational energy per unit mass 


E v - Ey/RT^ 


f 

G(T) 

H 

h 

kf 

k r 

k B 

K e 

m 

M 

n,k 

N 

[03 

p 


coupling factor defined by equation (17) 
term defined by equation (11) 
enthalpy per unit mass 
Planck constant 

forward reaction- rate constant, ce ^^ e ^ er j 

mole- second 

reverse reaction-rate constant, c ^-^ ieter ^ 

mole*- second 

Boltzmann constant 
equilibrium constant 
mass of an oxygen atom 
undisturbed free -stream Mach number 
number of increments in mesh system 

number of vibrational energy levels included in dissociation energy 
concentration of O2, mole/ centimeter^ 

pressure, atmospheres (1 atmosphere = 101.3 kilonewtons per meter2) 
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Qo’Qo 2 

R 

R 

t 

T 


electronic partition functions for O and 0 2 
gas constant for 0 2 , R = R/W 
universal gas constant 
time, seconds 

translational temperature, °K 


T = T/T. 


T 


m 


T 


v 


u,v 


W 


x,y 

a 

e 

®d’®R’®v 


P 

T 

Subscripts: 

e 


temperature defined by equation (12), °K 
vibrational temperature, °K 

velocity components parallel and normal to the wedge surface, respectively 
molecular weight of 0 2 

coordinates along and normal to wedge surface 
atom mass fraction 
wedge half-angle 

characteristic dissociation, rotational, and vibrational temperatures, 
respectively, °K 

mass density 

vibrational relaxation time, seconds 
local equilibrium conditions 


f 


forward reaction 



meas 


measured 


r reverse reaction 

vib eq vibrational equilibrium 

°° undisturbed free-stream condition 

Abbreviations: 

CVD coupled vibration-dissociation model of reference 1 

CVDV coupled vibration-dissociation-vibration model of reference 2 

ANALYSIS 

The gas model used here is a diatomic gas subject to harmonic vibrations, dissoci- 
ation, and recombination. Translation and molecular rotation modes are assumed to be 
fully excited, whereas the electronic internal energy and ionization are neglected. Since 
the gas model used for the wedge flows herein is patterned after that used for the normal 
shock wave computations in reference 2, the same dissociation reactions and vibrational 
and dissociational rates as used in reference 2 are used here. For more detail con- 
cerning the gas model, see references 1, 2, 3, 12, and 13. 

Figure 1 illustrates the flow geometry and body-oriented coordinate system for a 
wedge. The basic flow and rate equations for the steady, inviscid flow over a wedge are 
written in conservation form as (ref. 13): 


Continuity: 



8(pu) f 9(pv) _ Q 
8x 8y 

(1) 

x-momentum: 

8(p + pu 2 ) ^ 9(puv) _ Q 
8x 8y 

(2) 

y- momentum: 

8(puv) ] 8(p +pv 2 ) _ 0 
8x 8y 

(3) 

Vibrational energy rate: 

8(puE v ) 9(pv E v ) dE y 

8x 1 8y ^ dt 

(4) 
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Rate of dissociation: 


3(pua) 3 {pvoi) _ p dot 

3x 3y F dt 

(5) 

Energy: 


H + ^ * — = Constant 
2 

(6) 

The enthalpy and equation of state are 


H = ( 7 2 3Q ) RT + ^ ” «) E v + a ^®d 

(7) 

p = pRT(l + a) 

(8) 



For oxygen the vibrational rate is given by the CVDV model of Treanor and 
Marrone (ref. 2) as 


£ ■ - [ E ( T ' T v) - * C°< T > - at , r 

e(t,T v ) R0v |^exp^0 v /T m ) - 1 exp(N0 v /T m ) - 1 
G(T) = RQ v (^-zJj 


where 


(9) 


(10) 


(ID 



( 12 ) 


For oxygen, 


rp 


1.6 x 10" 9 exp f 10 *’^ 
>pl/3 


atm-sec 


N = 26 

0 V = 2274° K 
0 d = 59 380° K 


with the use of a linear harmonic oscillator, E v and E v>e can be written as (ref. 1): 

H®V J Tn ^®v 


E v = 


V " ©v/T v 


and E 


- 1 


v,e 


- 1 
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The vibrational energy rate for the CVD model is given by the simpler expression 


dE v _ E v e _ E v 
dt r 


( 13 ) 


The following reactions, which are also the reactions used in reference 2, are used for the 
dissociation of oxygen: 


Therefore, 


k f 

0 2 + X ^ 20 + X 


x = o, o 2 



^ w 

4p^o;3(l + q) 
W 2 (l - a) 


and the net dissociation rate for oxygen is 


da = i- p(l + oj) 
dt f W 


(1 - o) - 


4po;2 k r 

W kj 


(14) 


This system of equations is sufficient to calculate the flow-field properties, once 
appropriate expressions for the reaction-rate constants kf and k r have been deter- 
mined. Matthews (ref. 7) gives the measured dissociation- rate constant for oxygen as 


Mmeas - t 1 - 1 x 10 25 )t- 2-5 exp(- 


cm° 

mole-sec 


(15) 


Although the temperature range for the experiments of Matthews was 3000° K to 5000° K, 
equation (15) is used for much higher temperatures in the present applications. When the 
vibrational energy is in equilibrium with the translational mode, the forward and reverse 
rate constants are related through the equilibrium constant 

Kg = ( k f)vib eq (16) 

k r 


Since recombination for the flow over a wedge is negligible until the vibrational energy is 
close to equilibrium, the effect of vibrational nonequilibrium on k r is neglected. The 
equilibrium constant for the dissociation of oxygen reaction is (ref. 12) 

K = 4m Hk B T)3/2 Or/ c -0 v / T )Qq 2 e -0 d /T 
6 W W T V Qq 2 
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where 0 r = 2.059° K. As indicated in appendix B of reference 14, the electronic parti- 
tion functions for oxygen vary less than 5 percent for 1500° K < T < 8000° K. There- 
fore, average values over this temperature range are Qq^ =3.3 and Qq = 8.8. The 
coupling effect of vibrational nonequilibrium on dissociation for the CVD and CVDV 
models is given by reference 1 as 


k f _ 1 

1 - exp(-NO v /T m ) 

exp(0 v /T v ) - 1 

( k f)vib eq N 

exp(© v /T m ) - 1 

exp(© v /T) - 1 


(17) 


Equations (15), (16), and (17) can be used to compute kj and k r , provided a rela- 
tionship between kj and (kf) meas or e q and (kf) me as is known - A survey 

of the literature indicates that several investigations (such as refs. 15 to 17) assumed the 
dissociation-rate constant for oxygen was measured with the vibrational energy in equi- 
librium. In his shock-tube experiments with oxygen, Matthews (ref. 7) observed a small 
relaxation region aft of the shock wave, but upstream of the reaction region. He therefore 
associated this small relaxation region with the vibrational relaxation and concluded that 
the vibrational energy was in equilibrium in the reaction zone where the experimental 
values of the dissociation -rate constants were obtained. However, in a more recent 
report Treanor (ref. 10) explains this phenomenon very differently. Treanor indicates 
that aft of the shock wave the vibrational temperature rises rapidly from room tempera- 
ture then levels off at some higher temperature which is appreciably less than the trans- 
lational temperature. (See fig. 2.) It is in this "quasi-equilibrium" zone (where T v has 
leveled off) that the experimental values of dissociation- rate constants have been meas- 
ured. Although T v < T in the quasi-equilibrium zone, the magnitudes of the first two 
terms on the right-hand side of equation (9) are nearly equal and the third term (recom- 

bination) is negligible, hence ~ 0. Physically, the rate of increase of vibrational 

dt 

energy due to collisions is balanced by the loss of vibrational energy due to dissociation. 
Treanor further says, "thus during the entire time when the dissociation rate is being 
measured, the vibrational mode is out of equilibrium with the translational temperature, 
and the coupling of vibration and dissociation causes the observed rate to be less than the 
rate that would be obtained with local vibrational equilibrium." Although the quasi- 
equilibrium zone depicted in figure 2 is for 0.5 percent C >2 in 99.5 percent Ar, the normal 
shock computations in reference 2 indicate a similar quasi-equilibrium zone exists for 
undiluted C> 2 , which Matthews used in his experiments (ref. 7). 

In reference 10 Treanor also gives an iterative method to obtain (kf) v pj eq from 
(kf)meas* However, it is cumbersome to apply this technique to two-dimensional flow- 
field computations since it would require iterations at every mesh point. To illustrate 
the possible deviations in the flow-field properties for the nonequilibrium flow of oxygen 
over a wedge (using the CVD and CVDV couplings), two interpretations of measured 
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dissociation- rate constants are considered. The first interpretation assumes 
N meas ~ (^f)vib eq 

which is the assumption made in references 15 to 17. The other interpretation assumes 
that the measured rate constant already contains the vibrational nonequilibrium effects 
and thus sets 

(k-f) meas = ^f 

It is believed that the computed flow properties for the CVD and CVDV coupling models 
using any other method of obtaining kj from meas should lie between these two 
limiting conditions. 


Method for Numerical Computations 

Let the flow field be divided into an orthogonal mesh system as shown in the fol- 
lowing sketch: 



Then the coordinates of a general mesh point n,k are x = x Q + n Ax and 
y = k Ay. Approximate flow properties are first determined along the initial data line 
(n = 0) by using the exact wedge-tip gradients given in appendix C of reference 13. The 
numerical technique for computing the flow-field solution at increasing values of n is to 
replace the governing partial differential equations (eqs. (1) to (5)) with finite-difference 
equations. 

The system of partial differential equations given by equations (1) to (5) may be 
written as the one vector equation: 

^ + d m = 0 (m = 1, 2, 3, 4, 5) (18) 

The difference scheme, suggested by Lax (ref. 18), replaces the partial derivative 

( 3arr\ 

) k a * mes h P 0 ^ with the modified forward difference quotient 
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and replaces 



( a m)n+l,k " 2 [( a ni)n,k+l + (^jnjk-l] 
Ax 

with the symmetric difference quotient 

( b m)n,k+l ~ ( b m)n,k-l 
2 Ay 


Since only the mesh points n,k+l and n,k-l are involved at n for both of these dif- 
ference quotients, the term (d m ) n ^ is replaced by the average value 

\ [( d m)n,k+l + ( d m)n,k-l| 


With these substitutions the system of partial differential equations given by equation (18) 
is replaced with 


( a m)n+l,k 


2 [( a m)n,k+l + ( a m)n,k-l| “ 2 Ay[( br ^ n > k+1 " ( b m)n,k-lj 
~ ~2 [( d m)n,k+l + ( d m)n,k-l) ( m = ^ 


(19) 


In this manner the flow-field properties at n+l,k may be determined from those at 
n,k+l and n,k-l. Therefore, starting with the approximate flow properties along the 
initial data line (n = 0), the flow properties can be calculated at all points of the 
staggered mesh n = 1, 2, 3, . . . and k = 1, 3, 5, 7, . . . when n is odd; and 
k = 0, 2,4,6,. . . when n is even. (See the sketch.) 


Equation (19) cannot be used to compute the properties on the wedge surface (k = 0) 
unless imaginary properties are known at k = -1, which lies inside the wedge. It was 
found in reference 13 that the following difference scheme gave more accurate results 
along the surface than determining imaginary properties at k = -1 by an extrapolation 
or reflection technique and by using equation (19). Therefore, the finite-difference 
scheme for body points makes the following replacements: 



( a m)n+l,0 ~ ( a m)n-l,0 
2 Ax 




+ ( b m)n-l,oJ 


( d m)n,0 “* 


1 

2 



+ 
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Then the system of partial differential equations given by equation (18) becomes 


( a m + Ax d m)n+l,0 ( a m 


Ax d m)n-l,0 ” 


2 Ax 
Ay 




( 20 ) 


This equation can be simplified by noting that (b^^j q = (b ir J n _j g = 0 for m = 1, 2, 

4, and 5; and the m = 3 vector is not required on the surface because the boundary con- 
dition requires v = 0 at k = 0. Thus, the difference equation for the body points finally 
becomes 


( a m + Ax d m)n + l,0 = ( a m ' Ax d m)n-l,0 " ^ ( b m)n,l < m = 2 ’ 4 ’ 5 > (21) 

As a consequence of Lax's difference scheme an artificial viscosity is implicitly 
introduced and allows the computations to proceed downstream of the initial data line as 
if no shock wave were present at all. The shock wave appears in the solution, however, 
smeared over several mesh spaces while accurately giving the proper jump conditions 
across the shock. For more details concerning the numerical technique, consult refer- 
ence 13. 

Since the system of equations being solved is hyperbolic, the ratio of the mesh 
spacing in the x-direction to the mesh spacing in the y-direction is limited by the hyper- 
bolic stability criterion. The free stream as well as the shock layer must be considered 
in determining the stability criterion because the computations carry through the shock 
wave. For all the examples given here, the free-stream stability criterion was more 
stringent than that of the shock layer since the free-stream velocity is inclined at the 
wedge half-angle 0 to the coordinate system. 

RESULTS AND DISCUSSION 

It is well known that for the supersonic, inviscid flow of a frozen gas past a wedge, 
the attached shock wave is straight and the fluid properties are constant between the body 
surface and the wave. The same qualitative results also hold for a gas in equilibrium; 
however, in this case, the shock wave is closer to the body and the constant fluid proper- 
ties are different from those of frozen flow. In particular, the equilibrium pressure and 
temperature are less than the corresponding properties for frozen flow. 

When the conditions are such that nonequilibrium effects must be considered, the 
flow field is complicated. At the tip of the wedge the flow is frozen since it has just 
passed through the infinitesimally thin shock wave and has not had time to relax. The 
translational and rotational energy modes reach equilibrium almost immediately; 
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however, the vibrational and dissociational energies require a finite time to reach equilib- 
rium. For an undissociated free stream, the vibrational energy and the degree of disso- 
ciation increase along the surface of the wedge. These conditions in turn decrease the 
temperature and pressure but increase the velocity and density. The shock wave is 
inclined at the angle for frozen flow at the tip, but bends downward and approaches the 
equilibrium angle far downstream. As the flow along the surface achieves an equilibrium 
condition downstream, the pressure is the same as that for the equilibrium wedge solu- 
tion; however, the same is not true for the other flow variables. Downstream of the nose, 
three regions are observed in the profiles of properties between the body and the shock 
wave. Near the surface there exists an entropy layer since the streamlines passed 
through the steeper shock wave near the tip of the wedge. In the vicinity of the wave, the 
flow properties are close to being frozen since they have just passed through the shock 
wave and have not had time to relax. Between these two regions is a region which is 
essentially in equilibrium. This portion of the flow field has passed through the shock 
wave further upstream and thus has had time to relax to the equilibrium conditions. 

The cases given below were computed to determine the differences between the two 
limiting conditions of (kf) mea s = ( k ^) v ib eq and ( k f) me as = k f for the nonequilibrium 
flow of oxygen over a wedge. In each case the shock-wave angle at the tip of the wedge 
was chosen to be 59°, which is close to the shock detachment angle for the Mach numbers 
considered. The free-stream properties were chosen so that the shock strengths at the 
wedge tip would be the same as those for the normal shock solutions computed in refer- 
ence 2 for oxygen. 

Case I: 

M = 9.45099 0 = 42.379° 

Case II: 

M= 27.17159 0 = 43.361° 

For each case, 

T = 249° K p = 2.23 X 10" 4 atm 

= ® a °o = ® 

Figures 3 and 4 show the variation of T, E v , and a along the wedge surface for 
case I. Figure 4 illustrates the lag in the dissociational rate for (kf) meas = (kf) v ib e q 
as compared with (k^) meas = k f along the wedge surface. However, since the degree of 
dissociation is small, there is very little difference between the two limiting conditions 
for the variation of T and E v along the surface. Also, there is no noticeable 
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difference between the CVD and CVDV models for this case. Although T and E y are 
close to equilibrium at x = 80 cm, a is still increasing and additional computations 
showed that a had not achieved an equilibrium condition as far downstream as 
x = 4750 cm. 

The variation of T, E v , and a along the wedge surface for case II is shown in 
figures 5, 6, and 7, respectively. For this case there are significant differences between 
the two limiting conditions as well as between the CVD and CVDV models. Figure 6 
shows that the vibrational energy for the CVDV models levels off for 2 cm ~ x ~ 20 cm 
and then increases to an equilibrium value near x = 100 cm. This zone is similar to the 
quasi- equilibrium region noted in figure 2 for the flow aft of a shock wave. Figure 7 also 
depicts a lag in the dissociation rate near the tip of the wedge for (kf) meas = (kf) v ib e q 
as was found in case I. Figures 8, 9, and 10 show the profiles of T, E v , and a normal 
to the wedge surface at x = 3 cm. These figures illustrate that the shock-layer thickness 
is largest for the (kf) me as = (kf)vib eq condition with the CVDV model. 

The condition (kf) me as = ^f * or CVD model was found to yield the same results 
as computations with no coupling between the vibrations and dissociation (f = 1). This 
unusual circumstance is attributed to the fact that the coupling factor f given by equa- 
tion (17) only affects k r (for the (kf) meas = kf condition), and for the flow over a 
wedge recombination is only significant in the regions where the vibrational energy is 
near equilibrium (hence f « 1). 


CONCLUDING REMARKS 


To illustrate two different interpretations of measured dissociation- rate constants, 
the nonequilibrium flow of oxygen over a wedge was computed with the CVD and CVDV 
coupling models for the two limiting conditions (kf) me as = (^f)vib eq an< * (kf) m eas = ^f 
(where kf is forward reaction-rate constant and the subscripts meas and vib eq 
refer to measured and vibrational equilibrium, respectively). These limiting conditions 
should provide bounds for flow properties computed using any other method of obtaining 
kf from (kf) meas with the previously mentioned coupling models. 

For the lower Mach numbers, where dissociation is small, there is very little dif- 
ference between the flow properties for the two limiting conditions, with the exception of 
the atom mass fraction. There is a noticeable lag of a (atom mass fraction) with 
(kf) meas = (^f) vib eq behind « with (kf) mea s = kfj however, since the degree of disso- 
ciation is small, these effects are insignificant as far as the other flow-field properties 
are concerned. 

At the higher Mach numbers, where dissociation is significant, appreciable differ- 
ences are found between the two limiting conditions for all the flow-field properties. 
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Until a better understanding of the measured reaction rates is obtained, calculations 
should be made for the two limiting conditions to indicate the possible range in magnitude 
of the flow variables. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 20, 1966, 

129-01-08-18-23. 
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Figure 9.- Vibrational energy profile normal to wedge surface at x = 3 cm for case II. M = 27.17; 0= 43.36 




y, cm 


Figure 10.- Atom mass fraction profile normal to wedge surface at x = 3 cm for case II. M = 27.17; 6 - 43.36°. 


NASA- Langley, 1967 12 L“5269 



"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 
importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Scientific and technical information generated 

under a NASA contract or grant and considered an important contribution to 
existing knowledge. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information derived from or of value to NASA 
activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
Utilization Reports and Notes, and Technology Surveys. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



